AD-A123  237  NUMERICAL  STUDV  OF  PHASE  CONJUGATION  IN  STIMULATED  -  1/1 

BRILLOUIN  SCATTERING  FROM  AN  OPTICAL  HAVEGUIDE(U)  NAVAL 
RESEARCH  LAB  WASHINGTON  DC  R  H  LEHMBERG  28  DEC  82 
UNCLASSIFIED  NRL-MR-4985  F/G  28/6  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDAR0S-1B63-A 


NRL  Memorandum  Report  498S 


Numerical  Study  of  Phase  Conjugation 
in  Stimulated  Brillouin  Scattering 
from  an  Optical  Waveguide 

R  H.  Lehmberg 


MS*’®- 
Isi 


•".TO1 


' 


Laser  Plasma  Branch 
.  Plasma  Physics  Division 


December  30,  1982 


NAVAL  RESEARCH  LABORATORY 
Wuhlastoa.  D.C. 


rfT\G 


Approved  for  public  release;  distribuiion  unlimited 


SECURITY  CLASSIFICATION  of  This  page  I’R'N.n  D.l.  Entorod, 


REPORT  DOCUMENTATION  PAGE 


I.  report  number 


.Mr 


NRL  Memorandum  Report  4985 


4.  TITLE  Ca rd  Submit) 

NUMERICAL  STUDY  OF  PHASE  CONJUGATION  IN 
STIMULATED  BRILLOUIN  SCATTERING  FROM  AN 
OPTICAL  WAVEGUIDE 


7.  AuThORCj 

R.  H.  Lehmberg 


*.  PERFORMING  ORGANIZATION  NAME  and  address 

U.S.  Naval  Research  Laboratory 
Washington,  DC  20375 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

U.S.  Department  of  Energy 
Washington,  DC  20545 


read  instructions 

BEFORE  COMPLETING  FORM 


1  RECIPIENT’S  CATALOG  NUMBER 


S  TYPE  OF  REPORT  4  PERIOC  COVERED 


S.  PERFORMING  ORG  REPORT  NUMBER 


4  CONTRACT  OR  GRANT  numberoi 


to  PROGRAM  Element.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

DOE  AI08-79DP  40092(172) 
47-0859-0-2 


IZ  REPORT  DATE 


December  30,  1 082 

U  NUMBER  OF  PAGES 
26 


i 


MONITORING  AGENCY  NAME  4  ADD  RESSfU  dllltrtnt  from  Comrollln*  OllUo)  IS.  SECURITY  CLASS  to I  I  hit  roport) 

UNCLASSIFIED _ 

IS.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


14.  DISTRIBUTION  STATEMENT  fot  thlt  Roport, 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  the  ebetrect  entered  In  Block  20.  It  different  from  Report) 


19.  KEY  WORDS  (Continue  on  reveree  tide  II  neceeeery  end  Identity  by  block  number) 


Wavefront  conjugation  fidelity 
Aberration  compensation 
Stimulated  scattering 
Optical  waveguide 


Pump  depletion 


20.  ABSTRACT  (Continue  on  reveree  etde  II  n»e eeeery  end  Identity  by  block  number) 

Phase  conjugation  in  stimulated  Brillouin  scattering  (SBS)  from  an  optical  waveguide  is  studied 
numerically,  using  a  Cartesian  two  dimensional  propagation  code  (BOUNCE)  that  includes  pump 
depletion.  The  simulations  show  that  light  scattered  back  through  a  phase  aberrator  produces  a 
far  field  profile  closely  matching  that  of  the  incident  beam;  however,  the  near  field  intensity 
exhibits  large  and  rapid  spatial  inhomogeneities  across  the  entire  aberrator.  even  for  conjugation 
fidelities  H  as  high  as  98%^This  effect  can  be  only  partially  removed  by  spatial  filtering  techniques. 

The  dependence  of  H  upom'SBS  gain,  scattering  length  L,  and  the  average  intensity  T and 

_ (Continued) 


StCUAITV  ClASSiriCATion  pr  This  MCI  {•»*«*  Omtm  _ 

20.  ABSTRACT  (Continued) 

angular  divergence  dp  of  the  pump  beam  is  first  studied  in  the  low  reflectivity  regime  (i.e.,  ignor¬ 
ing  pump  depletion).  For  amplitude  gains  >el0,H  decreases  monotonically  with  the  factor 
1/9  D  due  to  a  small  scale  intensity  pulling  effect.  This  result  appears  to  be  independent  of 
whether  dD  arises  entirely  from  the  aberrator,  or  from  a  combination  of  the  abemtor  plus 
curvature  introduced  by  a  lens.  In  all  cases,  pump  depletion  is  found  to  enhance  the  fidelity  by 
inhibiting  the  small  scale  pulling  effect.  Under  appropriate  conditions,  this  can  actually  reverse 
the  intensity  dependence  of  H  seen  in  the  absence  of  pump  depletion. 


steumrr  classification  or  this  VAoefWfcM  om>  i*im« 


ii 


CONTENTS 


i  -  ■  ■  *•  .  . 

I.  INTRODUCTION  . I 

II.  THEORY  AND  CODE  DESCRIPTION  .  ^2 

III.  LOW  REFLECTIVITY  LIMIT  .  4 

A.  Aberrated  Plane  Wave  .  4 

B.  Aberrated  Diverging  Wave  .  7 

IV.  EFFECT  OF  PUMP  DEPLETION .  7 

V.  DISCUSSION  .  8 

ACKNOWLEDGMENTS  .  9 

REFERENCES  .  9 


iii 


NUMERICAL  STUDY  OF  PHASE  CONJUGATION 
IN  STIMULATED  BRILLOUIN  SCATTERING 
FROM  AN  OPTICAL  WAVEGUIDE 


I.  INTRODUCTION 

Optical  wavefront  conjugation  in  stimulated  backscatter  has  been  the  subject  of  extensive  experi¬ 
mental  and  theoretical  investigation  over  the  last  ten  years.'  Nevertheless,  there  remain  significant  gaps 
in  the  quantitative  understanding  of  this  phenomenon.  Analytic  theories2-10  generally  agree  that  it 
occurs  because  the  component  of  backscattered  light  conjugate  to  a  spatially  inhomogeneous  pump 
wave  will  grow  at  approximately  twice  the  rate  of  any  random  modes  that  are  orthogonal  to  it,  and 
should  presumably  dominate  at  sufficiently  high  gain.  They  disagree,  however,  on  the  question  of  why 
some  of  the  backscatter  remains  unconjugated,  even  in  the  high  gain  regime  where  the  conjugate  frac¬ 
tion  or  fidelity  H  should  approach  unity.  It  is  an  important  question  because  even  a  small  residue  of 
unconjugated  light  (e.g.  1  —  H  ~  1%)  can  cause  serious  beam  distortion  in  the  near  field  of  the  abera- 
tor.  Zel'dovich4  5  links  this  residue  to  small  scale  intensity_"pul!ing"  by_the  inhomogeneous  pump 
beam,  and  predicts  that  1  —  //should  increase  with  the  factor  7/0 5,  where  7 is  the  average  intensity  and 
0O  is  the  angular  divergence  of  the  pump  light.  However,  since  the  analysis  is  based  on  plane  wave 
expansions,  it  may  be  inaccurate  for  describing  waveguides.  Heilwarth6  uses  the  correct  waveguide 
modes,  but  ignores  the  terms  responsible  for  the  pulling  effect.  He  predicts  that  under  high  gain  condi¬ 
tions,  the  fidelity  should  depend  upon  the  distribution  of  pump  power  among  the  modes  (i.e.  upon  the 
detailed  spatial  characteristics  of  the  aberration),  but  not  upon  /.  With  the  exception  of  Ref.  (10). 
none  of  these  theories  includes  pump  depletion,  which  becomes  important  when  the  backscatter  is 
significantly  above  threshold.  The  main  result  in  Ref.  (10)  is  that  pump  depletion  will  drastically  limit 
the  allowable  gain,  and  hence  the  conjugation  fidelity,  if  the  initial  backward  propagating  noise  level  is 
too  high.  This  condition  is  not  applicable  to  the  operating  regime  of  interest  here. 

In  a  recent  paper,"  we  described  a  numerical  study  of  wavefront  conjugation  in  stimulated  Bril- 
louin  scattering  (SBS),  using  a  steady  state  2D  propagation  code  (BOUNCE)  that  includes  pump  deple¬ 
tion.  These  calculations  simulated  both  aberration  correction1 2-21  and  image  replication22-24  experi¬ 
ments,  but  they  were  restricted  to  the  case  of  a  beam  focused  into  a  medium  unbounded  in  the 
transverse  direction.  The  present  work  examines  SBS  in  a  multimode  optical  waveguide,  and  studies 
the  parameters  that  affect  the  wavefront  conjugation  fidelity.  The  key  results  are  summarized  below. 

(i)  The  waveguide  was  able  to  produce  significantly  higher  fidelities  ( H  <  98%)  than  the 
focused  configuration  (H  <  88%),  in  agreement  with  several  experimental  studies.15-19 
This  improvement  occurs  because  waveguides  can  avoid  the  deleterious  spatial  gain  nar¬ 
rowing  effects  that  invariably  accompany  Gaussian-like  pump  beam  profiles  near 
focus.®'9, 11  ■ 25 

(ii)  The  simulations  show  that  the  light  scattered  back  through  the  phase  aberrator  exhibits  a 
far  field  intensity  profile  closely  matching  that  of  the  unaberrated  incident  beam,  as  seen 
in  experiments,12-21  while  the  unconjugated  portion  appears  mainly  as  low  level  hash  dis¬ 
tributed  over  large  angles  outside  this  profile.  However,  the  near  field  intensity  exhibited 
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large  and  rapid  spatial  inhomogeneities  across  the  entire  aberrator  in  ail  cases,  even  when 
the  phase  remained  uniform  to  within  X/15,  and  the  fidelity  was  98%.  Similar  results 
were  found  in  the  previous  work.11  An  attempt  was  made  to  remove  the  inhomogeneities 
by  a  spatial  filtering  technique,  but  this  was  only  partially  successful. 

(iii)  In  the  absence  of  pump  depletion,  the  fidelity  was  found  to  increase  with  the  average 
pump  intensity  for  amplitude  gains  up  to  around  e10,  then  decrease  slowly  and  monotoni- 
cally  with  higher  intensity.  The  low  intensity  behavior  agrees  with  the  usual  mode 
theory,2-7  which  attributes  SBS  wavefront  conjugation  to  the  higher  gain  experienced  by 
the  conjugate  wave.  At  high  intensity,  the  beam  profiles  dearly  show  a  small  scale  pulling 
effect,  in  qualitative  agreement  with  Zel’dovich;4,5  however,  the  detailed  parameter  study 
presented  here  suggests  that  If  varies  as  9D/7,  rather  than  the  9 1/7  predicted  in  Refs.  (4) 
and  (5).  A  heuristic  argument  linking  this  to  the  discrete  nature  of  the  waveguide  modes 
is  presented  in  Sec.  V.  All  of  these  results  appear  to  be  independent  of  whether  9D  arises 
entirely  from  an  aberrator,  or  from  a  combination  of  the  aberrator  plus  phase  curvature 
introduced  by  a  lens.  This  explains  why  it  has  been  possible  to  achieve  high  conjugation 
fidelities  in  a  waveguide  with  only  a  moderately  aberrated  (or  even  unaberrated)  pump 
beam.16 

(iv)  In  all  of  the  cases  studied,  pump  depletion  significantly  enhanced  the  fidelity  of  the  wave- 
front  conjugation  by  inhibiting  the  small  scale  pulling  effect.  A  similar  result  was  found 
in  the  focused  configuration.11  If  one  chooses  conditions  such  that  pump  depletion  starts 
to  become  significant  at  gains  around  e10.  then  the  fidelity  continues  to  increase  with  the 
incident  intensity,  rather  than  decreasing  as  in  the  cases  described  above.  This  result  is  in 
good  qualitative  agreement  with  the  experiments  of  Mays  and  Lysiak.19 

Section  II  gives  a  brief  description  of  the  nonlinear  propagation  code,  and  defines  the  quantities 
that  are  calculated.  Section  III  describes  the  parameter  study  in  the  low  reflectivity  limit,  while  Sec.  IV 
looks  at  the  effects  of  pump  depletion.  These  results  are  summarized  and  discussed  in  Sec.  V. 

II.  THEORY  AND  CODE  DESCRIPTION 


Under  steady  state  conditions,  the  complex  pump  and  SBS  backscatter  amplitudes  EL  (x.  r)  and 
£$(x  z)  satisfy  the  parabolic  equations 


d  id2 
dz  2k  dx2 


El 


jd _ i_  d2  p 

dz  2k  3x2  5 


(la) 


(lb) 


in  a  two  dimensional  cartesian  geometry.  Here,  k  —  2nn/\  is  the  magnitude  of  the  propagation  vectors 
(assumed  equal  in  this  paper),  x  is  the  transverse  coordinate,  and  g  is  the  coupling  constant  of  the  Bril- 
louin  medium,  which  is  contained  within  a  region  zx  <  z  <  z2.  For  simplicity,  we  will  take  n  s  1, 
and  ignore  refraction  in  this  paper.  BOUNCE  solves  Eqs.  (la,b).  assuming  an  aberrated  pump  wave 
El(x,  z2)  is  incident  at  z2,  while  the  backscatter  grows  from  a  small  counterpropagating  noise  wave 
Es(x,  Z|)  introduced  at  zx.  (E.g.,  see  Figs.  2  and  8.) 

Since  £4  and  Es  are  specified  on  opposite  ends  of  the  medium,  and  neither  steady  state  solution  is 
known  within  that  medium  at  the  outset,  the  inclusion  of  pump  depletion  is  a  nontrivial  problem.  One 
must  know  Es(x,  z)  in  order  to  calculate  EL (x.  z),  and  vice  versa.  The  brute  force  solution  would  be 
to  simply  integrate  the  time  dependent  equations  over  many  optical  transits  through  the  medium  until 


steady  state  is  finally  attained.  The  approach  followed  here  is  an  iteration  procedure,  which  is 
significantly  faster.  Basically,  it  starts  with  an  undepleted  pump  wave  E[lHx,  z),  calculates  the 
resultant  backscatter  £5"  Or.  z)  throughout  the  medium,  then  alternately  recalculates  these  waves  until 
the  solutions  converge  to  self  consistency.  A  detailed  description  of  the  techniques  used  to  effect  this 
convergence  is  given  in  Ref.  (11). 

In  the  previous  work,  all  of  the  beam  profiles  were  assumed  to  be  contained  entirely  within  the 
transverse  region  -  W/ 2  <  .*  <  +  fV/2,  taking  care  to  ensure  that  empty  zones  were  maintained 
around  ±  W/2  to  simulate  propagation  in  an  unbounded  medium.  The  calculations  were  then  per¬ 
formed  by  a  split  operator  technique,26  which  expanded  the  field  amplitudes  using  the  conventional 
discrete  Fourier  transform  (DFT)  with  N  transverse  points: 

E(J)  -  jJ  E(K)  exp(i2irJK/N)  (2a) 

A>0 

xU)  -  (W/N)J-  W/2,  J-  0.1 . tf-l  (2b) 

This  expansion  has  periodic  boundary  conditions:  hence,  if  a  wave  were  allowed  to  propagate  to  one 
edge  (e.g.,  +  W/2)  it  would  re-enter  from  the  opposite  side  (e.g.,  —  W/2)  as  illustrated  in  Figs.  (la-c). 

In  this  paper,  we  have  simulated  propagation  in  an  optical  waveguide  of  width  W,  where  the  field 
amplitudes  satisfy  the  boundary  conditions 

Els(  W/2,  z)  -  Els(-  WI2,  z)  -  0  (3) 

in  all  cases;  hence,  waves  will  reflect  from  the  boundaries,  as  illustrated  in  Figs.  (Id-f).  This  can  be 
accomplished  by  replacing  the  usual  DFT  by  a  waveguide  mode  expansion 

EU)  -  y  E(K)sin(nJK/S).  J  -  0.  1  ....  N  -  I.  (4) 

*-i 

which  automatically  satisfies  (3).  In  the  simulations  presented  here,  we  followed  an  equivalent 
approach,  which  uses  the  conventional  DFT  with  2 N  points.  The  additional  .V  points 
N,  N  +  1 . 2.V  —  1  describe  the  propagation  of  an  external  "image"  field  that  satisfies 

E’(2N-  /)-  —  EU),  /-  1,2 . AT-1.  (5a) 

£'(AO  -  £(0)  -  0.  (5b) 

If  the  conventional  DFT  (with  2N  points)  is  applied  to  these  combined  waves,  it  will  populate  only  the 
desired  waveguide  modes  sin(nJK/N). 

The  accuracy  of  the  modified  propagation  code  was  checked  in  the  following  ways: 

(i)  Reruns  of  some  of  the  earlier  simulations,  in  which  the  beams  remained  isolated  from 
±  W/2,  were  found  to  give  identical  results. 

(ii)  To  test  free  propagation  within  the  waveguide,  the  numerical  results  were  compared  to  an 
analytic  calculation  (based  on  multiple  images)  for  the  case  where  an  unaberrated  cylindri¬ 
cal^  diverging  wave  was  introduced  at  the  entrance  plane.  The  agreement  was  excellent. 

(iii)  If  pump  depletion  is  ignored,  the  total  number  of  points  ordinarily  used  in  the  calculation 
(256  x  601)  can  be  doubled,  thus  allowing  a  finer  grid  in  either  .tor  z.  When  several  of 
the  simulations  were  rerun  with  these  finer  grids,  the  agreement  was  in  all  cases  better 
than  1%  for  the  calculated  values  of  H. 

The  code  evaluates  several  important  parameters  in  addition  to  the  usual  pump  and  backscatter 
beam  profiles.  These  include  the  power  gain  G  and  total  reflectivity  R: 
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G  s  Ps(z2)/Ps(z ,),  R  &  Ps(z2)/PL(z2),  (6a, b) 

where 

Pl.s (2)  =  /  lfi.sU  z)l2  <*r  (7) 

is  proportional  to  the  total  pump  or  backscatter  power,  assuming  unit  width  along  the  y  direction.  From 
Eqs.  (la,b)  and  (3),  one  can  readily  derive  the  Manley-Rowe  relation 

PL(z)  -  Ps(z )  —  constant  (8) 

for  <*s  —  <aL~  The  code  also  calculates  the  phase  conjugate  fidelity  defined  by  the  normalized  correla¬ 
tion  function 

Hiz)  S  PL(z)Ps(z)  1  S-wn  £*U  z)Es(x'  z)  (9) 

(Note  that  H  —  l  in  the  case  of  perfect  conjugation  Es  «  El).  According  to  Eqs.  (la,b)  and  (3),  H 
remains  constant  when  g  *  0;  thus  H  -  H (r2)  at  all  points  beyond  z2. 

III.  LOW  REFLECTIVITY  LIMIT 

This  section  will  deal  with  the  case  where  the  reflectivity  is  low  enough  that  pump  depletion  can 
be  ignored,  and  thus  the  pump  power  PL  will  be  constant  along  z 

A.  Aberrated  Plane  Wave 

In  order  to  obtain  a  direct  comparison  with  analytic  theories,  we  first  consider  the  "pure"  case 
where  the  angular  beam  divergence  ft0  arises  only  from  the  aberrator  (ftD  -  ft.,).  as  shown  in  Fig.  2.  A 
uniform  plane  wave  of  k  —  694  ^m  traverses  a  hard  aperture  of  width  W  —  2  mm  followed  by  aberra¬ 
tor  A ,  which  imposes  a  random  phase  modulation  and  concomitant  far  field  angular  divergence  ft ,  20 

mrad  (full  1/e  width),  as  shown  in  Fig.  3.  The  aberrated  beam  ( — 60  x  diffraction  limited)  then 
immediately  enters  the  optical  waveguide,  also  of  width  W. 

In  all  of  the  simulations  presented  here,  the  beam  was  allowed  to  propagate  20  cm  into  the 
waveguide  before  entering  the  scattering  medium  z2.  This  was  done  in  order  to  ensure  a  complete  mix¬ 
ing  of  radiation  from  different  parts  of  the  aberrator,  so  that  pronounced  intensity  (as  well  as  phase) 
inhomogeneities  were  present  within  the  entire  gain  length  L ;  again  for  purposes  of  comparison  with 
the  analytic  theories.  In  spite  of  these  small  scale  inhomogeneities,  the  intensity  will  remain  statistically 
uniform  across  the  waveguide.  This  is  in  contrast  to  the  focused  configuration,11  where  the  envelope  of 
the  intensity  profile  was  typically  Gaussian-like.  Test  runs,  in  which  the  free  propagation  distance  was 
eliminated,  yielded  similar  results,  but  there  was  more  statistical  spread  in  the  data  with  different 
members  of  the  aberrator  ensemble  (i.e.,  different  statistical  realizations). 

The  noise  source  at  zx  was  modelled  by  a  random  complex  field  Es(x,  ;,)  of  100  mrad  angular 
divergence,  giving  initial  pump-noise  correlations  H t  -  H(z\)  typically  in  the  range  0.001  to  0.01. 
However,  the  wavefront  conjugation  of  the  backscatter  was  found  to  be  generally  insensitive  to  the  sta¬ 
tistical  properties  of  the  noise  source  (including  even  a  coherent  source)  as  long  as  the  gain  was  high 
enough  and  H\  remained  small. 

Referring  now  to  the  lower  left  hand  corner  of  Fig.  2,  one  sees  that  as  Es(x,  ;)  grows  and  pro¬ 
pagates  toward  z2,  it  develops  an  intensity  distribution  with  small  scale  inhomogeneities  similar  to  those 
of  |£t,Gr,  z)|2,  giving  a  fidelity  H  —  93%.  From  z2>  the  backscatter  retraces  the  pump  beam  path 
through  the  aberrator  to  produce  the  near  and  far  field  intensity  profiles  shown  on  the  right  hand  side 
of  Fig.  2.  The  phase  compensation  is  evident  in  the  far  field,  where  the  backscatter  exhibits  a  strong 
central  spike  (of  width  2k/  W)  and  nearby  sidelobes  almost  identical  to  the  sinc(.v)  angular  spectrum  of 


the  unaberrated  incident  beam.  The  —7%  unconjugated  component  appears  mainly  as  low  level  back¬ 
ground  hash  extending  well  beyond  this  central  structure;  hence,  expression  (9)  gives 

l  -  ff  ^  background  hash 

total  backscatter  power  ‘ 

For  example,  if  the  integrals  in  expression  (9)  were  evaluated  only  over  the  central  spike  and  first  three 
sidelobes.  H  would  jump  to  99.3%.  Although  the  peak  levels  of  the  background  are  typically  less  than 
1%  of  the  on-axis  intensity,  they  are  still  significantly  larger  than  the  experimentally  observed 
values.15-14  This  stems  from  the  fact  that  the  hash  is  confined  to  only  one  transverse  dimension;  i.e..  it 
lacks  the  usual  geometrical  weighting  factor  that  would  require  lower  levels  in  the  analogous  two- 
dimensional  "halo"  in  order  to  satisfy  expression  (10). 

In  the  near  field  of  the  aberrator,  the  conjugated  and  unconjugated  components  become  mixed, 
resulting  in  serious  degradation  of  the  entire  profile.  This  behavior,  which  was  found  in  all  cases  both 
here  and  in  the  focused  configuration,11  has  gone  largely  unnoticed  in  the  aberration  correction  experi¬ 
ments,  since  these  normally  record  only  far  field  intensities.  It  will  be  discussed  in  greater  detail  in  Sec. 
V. 


Although  the  analytic  theories2-6  disagree  significantly  on  some  issues,  they  do  agree  on  the 
primary  mechanism  that  leads  to  wavefront  conjugation;  i.e.,  that  the  component  of  the  backscatter 
radiation  proportional  to  £|  (x  r)  grows  at  approximately  twice  the  rate  of  the  non-conjugated  modes. 
Since  the  conjugate  field  ££(x  z)  satisfies  the  equation 


dz 


2  k  dx2 


El(x.  z)  -  0 


in  the  absence  of  pump  depletion,  the  lowest  order  solution  of  Eq.  (lb)  can 
approximation1-6- 12 


Es(x,  z)  —  C  E*l(x,  z)  exp 


~crgi(z  -  Z|) 


+  £vc(x  z)  exp 


lg 


where  C  is  a  constant,  /  -  PJW is  the  average  pump  intensity,  and  £Vc(x 
Eq.  (11)  that  are  orthogonal  to  E*L(x.  z): 


w/2 

S-wn  El (x  £vc('v'  :)dx  "  °- 


(11) 

be  written  in  the  mode 

7(z  -  r,)J.  <12) 

is  a  sum  over  modes  of 

(13) 


The  theory  then  predicts  that  <r  —  2;  hence,  the  first  term  of  (12)  will  eventually  predominate  if 
g7(z  -  Zi)  becomes  large  enough,  even  though  most  of  the  Stokes  energy  at  Z|  resides  in  the  second 
term.  Physically,  the  condition  <r  >  1  is  linked  to  the  multimode  intensity  structure  of  the  pump 
beam.12 


One  can  relate  this  theory  to  the  numerical  calculations  by  substituting  (12)  into  expression  (9), 


and  using  (13)  to  obtain 

Hy  PS(Zy)  -  |C|2Pt 

(14a) 

///»s(z2)-  |C|2Pt  expWL). 

(14b) 

where  H\  -  H(zx),  H  -  H(z2),  and  L-  z2- 

Z|.  Definition  (6a)  then  yields 
ln(GH/Hy) 

(15) 

<T  ” 

glL 

The  numerical  calculations  performed  here  found  1.95  <  <r  ^  2.4,  in  approximate  agreement  with 
mode  theory  and  recent  measurements.20  Surprisingly,  the  largest  values  of  a-  occurred  at  low  fidelities. 
This  is  because  the  large  o-  and  reduction  of  H  are  caused  by  the  same  intensity  pulling  effect,  which 
will  be  discussed  in  greater  detail  below. 
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A  more  revealing  measure  of  the  limitations  of  the  mode  approximation  is  provided  by  the  ratio  p 
of  output  power  in  the  conjugated  and  unconjugated  components  of  (12): 

_  |C|2  PL  cxp(<rglL) 

p  st  - = - ,  (16) 

P»c  exp (glL ) 

where 

W/2 

Psc  =  Lwn  {£ncU  =‘)l2  dx-  (l7) 

To  rewrite  (16)  in  terms  of  the  calculated  quantities,  we  first  substitute  (12)  into  (7)  and  use  (14a)  to 
obtain 

P\c  -  Ps<*i>  ~  -  (1  -  //,)Fs<z,)  =  Psizt):  (18) 

then  combine  this  result  with  (14b)  and  (16): 

p  -  GH/expigJL).  (19) 

Figure  4  shows  a  plot  of  Inp  and  the  fidelity  H  as  functions  of  the  average  gain  increment  glL, 
using  tjie  same  20  mrad  aberrator  and  7S  cm  length  used  in  Fig.  2.  At  low  gains,  the  fidelity  improves 
with  g/L,  as  one  would  expect,  since  p  is  also  increasing.  With  further  increase  of  git.  however,  H 
goes  through  a  broad  maximum,  then  slowly  declines  instead  of  approaching  unity  for  Inp  »  1,  as  Eq. 

(12)  would  predict.  Similar  behavior  was  observed  with  other  aberrators  and  waveguide  lengths,  with 
the  maximum  generally  occurring  in  the  range  9  <  GlmML  <  11.  (As  expected,  the  larger  values  of 
7 usually  occurred  in  cases  where  the  initial  correlation  H\  was  smaller.)  This  behavior  occurs 
because  the  inhomogeneous  pump  beam  tends  to  "pull"  a  disproportionate  amount  of  the  backscattered 
power  into  the  most  intense  spikes,4  3  as  seen  in  the  ;2  Plots  in  Figs.  2  and  5.  It  is  therefore  a  small 
scale  analog  of  the  spatial  gain  narrowing  effect  seen  in  the  focused  configuration.11  Since  the  pulling 
tendency  arises  from  the  exponential  nature  of  the  gain  term  in  Eq.  (lb),  and  is  counteracted  only  by 
the  diffraction  term,  it  is  expected  to  become  more  evident  as  the  gain  e-folding  length/diffraction 
length  ratio  kO^/gl  —  kLOo/g/L  decreases.  This  is  illustrated  qualitatively  by  comparing  Figs.  2  and  5, 
which  use  the  same  20  mrad  aberrator  and  approximately  the  same  gain  exponent  glL .  but  differ  by  a 
factor  of  3.75  in  length. 

Zel'dovich  and  coworkers4, 5  have  attempted  to  treat  the  pulling  effect  quantitatively,  using  a  per¬ 
turbation  analysis  with  the  mode  approximation  as  zeroth  order.  They  found  that  H  depends  only  on 
the  single  parameter  key g7,  approaching  unity  asymptotically  for  key gl  »  1.  Here,  the  2D  version 
of  this  theory5  is  compared  to  the  numerical  calculations  from  BOUNCE,  using  different  values  of  L 
and  9D.  For  each  value  of  L  and  90,  calculations  were  performed  on  five  different  members  of  the 
aberrator  ensemble,  in  order  to  ensure  an  adequate  statistical  sampling.  Following  Zel’dovich,  we 
always  kept  the  total  gain  glL  high  enough  to  discriminate  against  the  non-conjugated  components  in 
the  mode  approximation.  In  practice,  this  meant  keeping  In(GH)  -  glL  >  5,  which  resulted  in  choos¬ 
ing  gains  somewhat  higher  than  those  giving  the  maximum  fidelity.  (Recall  Fig.  4.)  Figure  6  compares 
the  results  from  Ref.  5  with  numerical  runs  made  for  90  -  20  mrad  (open  diamonds)  using  L  «  20, 

40,  and  75  cm,  and  9D  —  10  mrad  (solid  diamonds)  using  L  -  40,  75,  200,  and  300  cm.  The  disagree¬ 
ment  between  the  theories  is  apparent  in  spite  of  the  spread  in  the  numerical  data.  In  Fig,  7,  these 
numerical  results  are  replotted  as  a  function  of  k0D/gl along  with  additional  data  for  9D  «  5  mrad  * 

(open  squares)  usingL  -  200  and  400  cm.  We  then  see  that  H  appears  to  be  a  function  of  k90/gl ’  r 

instead  of  the  koygJ predicted  in  Refs.  (2-5),  at  least  for  the  range  of  parameters  studied  here.  The 
probable  reason  for  the  disagreement  is  that  those  treatments  are  based  upon  an  expansion  in  plane 
waves  rather  than  waveguide  modes,6  a  point  that  will  be  explored  in  greater  detail  in  Sec.  V. 
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B.  Aberrated  Diverging  Wave 

So  far,  the  results  have  been  restricted  to  the  "pure"  case  where  the  beam  divergence  arises 
entirely  from  the  aberrator;  i.e.,  where  9D  -  9A.  However,  in  most  phase  compensation  experiments 
that  use  an  optical  waveguide1213, 15-20  the  aberrator  is  imaged  onto  or  near  the  entrance  face  of  the 
guide,  as  illustrated  in  Fig.  8.  The  light  emerges  from  this  image  with  a  spherical  (or  in  this  case 
cylindrical)  wavefront  curvature  of  radius  f/M,  where  /  is  the  focal  length  of  the  lens,  M  5  WJ  W  is 
the  magnification,  and  W4  is  the  aperature  at  the  aberrator.  This  corresponds  to  a  maximum  ray  diver¬ 
gence  angle 

9F-  WMlf-  Wjf  (20) 

in  addition  to  the  random  divergence  0'4  -  M9A  due  to  the  aberrator.  After  the  light  propagates  within 
the  guide  a  distance  of  order  W/{9F  4-  0'4)  beyond  the  image,  these  components  should  mix  to  produce 
an  effective  random  beam  divergence  9D  satisfying  the  condition 

max(9F,  9A )  <  9D  <  0F  +  9'4. 

Several  numerical  runs  were  made  with  the  configuration  shown  in  Fig.  8,  using  9F  * 
and  9'a  —  5  mrad  (•— 15  x  diffraction  limit).  As  in  the  "pure"  case,  a  20  cm  spacing  was  left 
and  the  aberrator  image  at  the  waveguide  entrance.  Within  the  context  of  the  mode  approxinr 
(IS)  again  yielded  <r  values  in  the  range  1.9S-2.4,  whUe  the  fidelity  showed  a  gain  dependet 
to  the  one  in  Fig.  4,  with  a  broad  maximum  around  g/L  =  10—12  and  Inp  >  5  just  beyonr 
imum. 

Following  the  same  procedure  as  before,  we  again  plotted  the  calculated  fidelity  vs  k9D/gl,  where 
9d  was  arbitrarily  chosen  to  be  the  vector  sum  9D  S  (9F  +  0^2),/2  =  20.6  mrad.  These  results  (shown 
by  the_asterisks)  are  compared  to  the  "pure"  cases  in  Fig.  9.  The  generally  good  agreement,  at  least  for 
k9D/gI  >  5  x  10J,  tends  to  confirm  the  mixing  hypothesis  and  the  important  contribution  of  9F  to  the 
effective  random  divergence  &D  within  the  gain  medium.  The  configuration  modelled  here  is  typical  of 
most  optical  waveguide  SBS  experiments,  where  9F  »  9'4>  and  thus  9D  =  0F  is  only  weakly  dependent 
upon  94.  This  explains  why  the  fidelity  has  been  found  generally  insensitive  to  the  aberrator  diver¬ 
gence  in  these  experiments.16 

If  the  aberrator  is  removed  entirely  fromjhe  cylinrically  diverging  beam,  //follows  a  curve  similar 
to  Fig.  9,  but  it  rises  more  rapidly  with  k90/gl\  e.g.,  the  emulations  gave  H  *  20%  at  k9D/gl  —  2900, 
86%  at  5900  and  98%  at  104.  There  is  no  reason  to  expect  the  usual  analytic  theories  to  apply  in  this 
case,  however,  because  there  is  no  longer  a  random  mix  of  different  portions  of  the  incident  beam. 

IV.  EFFECT  OF  PUMP  DEPLETION 

This  section  deals  with  the  case  where  the  reflectivity  R  is  sufficient  that  pump  depletion  can 
become  an  important  factor.  We  will  examine  only  the  configuration  shown  in  Fig.  8,  since  this  relates 
more  closely  to  actual  experiments. 

Figure  10  shows  the  reflectivity  and  fidelity  vs  the  average  incident  intensity  /(ri)  for  L  ™  75  cm. 
9f  *■  20  mrad,  and  one  member  of  the  ensemble  of  aberrators  with  9'4  -  M9 4  ■*  5  mrad.  _The  dashed 
line  shows  the  fidelity  that  would  be  observed  if  pump  depletion  were  ignored:  i.e.,  if  /(;)  -  /(::) 
were  valid  throughout  the  length  of  the  gain  medium.  It  closely  resembles  the  corresponding  curve  in 
Fig.  4,  with  a  maximum  around  gI(z2)L  —  11.  With  pump  depletion,  however,  //continues  to  grow 
with  increasing  intensity  and  reflectivity,  in  agreement  with  the  experimental  results  of  Mays  and 


(21) 

mrad* 
een  z2 
->n,  Eq. 

imilar 
-■  ax- 


*  Allhough  calculations  with  larger  divergence  angles  would  be  of  interest,  the  additional  memory  storage  required  to  ensure  accu¬ 
racy  is  somewhat  beyond  the  present  capabilities  of  the  NRL  computing  system. 


Lysiak.19  This  occurs  because  the  large  backscatter  tends  to  counteract  the  pulling  effect  by  selectively 
depleting  the  more  intense  spikes  in  the  pump  beam  at  z  <  r2.  An  analogous  effect  was  also  seen  in 
the  focused  configuration.11  It  is  worth  emphasizing  here  that  if  pump  depletion  is  present,  the 
effective  average  gain  increment  is 

y  =  -f  P  7<z)rfz  =  (l  -  R)g7(z2).  (22) 

which  can  be  significantly  smaller  than  gi(z2)\  e.g.  at  giiz2)L  =*  22.5  (corresponding  to  R  =  41%),  y 
was  calculated  directly  to  be  13.5,  in  good  agreement  with  (1  -  R)g/(z2)L  »  13.3. 

No  attempt  was  made  to  plot  H  vs  kOD/gI(z2)  for  the  finite  depletion  case.  (Such  calculations 
would  be  costly  because  each  point  requires  10-20  iterations.)  For  R  <  10%,  however,  the  results 
shown  in  Fig.  9  should  remain  a  reasonable  (if  slightly  pessimistic)  approximation,  while  at  higher 
reflectivities,  the  fidelity  should  at  least  maintain  a  similar  functional  dependence  upon  0D. 

Figure  11  shows  the  beam  profiles  corresponding  to  the  41%  reflectivity  point  of  Fig.  10.  In  spite 
of  a  high  conjugation  fidelity  (H  =  97%),  which  is  evident  in  both  the  near  field  phase  and  far  field 
intensity  profiles,  the  near  field  intensity  profile  continues  to  exhibit  large  random  inhomogeneities.  An 
attempt  was  made  to  eliminate  this  problem  by  spatial  filtering,  using  far  field  apertures  of  lO0t  and 
40t,  where  0t  —  2 \/WA  is  the  width  of  the  central  far  field  spike.  As  seen  in  Fig.  12,  this  approach 
was  only  partly  successful:  thus,  even  the  small  unconjugated  components  in  the  near-axial  sidelobes 
(<  1%)  can  significantly  distort  the  near  field. 

V.  DISCUSSION 

The  numerical  calculations  described  in  Sec.  III  have  explored  the  effects  of  the  average  gain  g7L. 
gain  increment  gl  and  angular  divergence  0O  upon  the  conjugate  fidelity  H  in  the  limit  of  negligible 
pump  depletion.  In  particular,  they  have  studied  the  limitations  of  the  mode  approximation  (and  con¬ 
sequent  loss  of  fidelity)  arising  from  the  small  scale  pulling  effect.  The  resulting  dependence  of  H 
upon  0O,  rather  than  the  05  predicted  by  analytic  theory,4  5  appears  to  be  related  to  the  discrete  nature 
of  the  waveguide  modes.  Although  a  detailed  treatment  of  the  problem  is  beyond  the  scope  of  this 
paper,  the  result  can  be  at  least  qualitatively  understood  by  the  following  heuristi  argument. 
Zel’dovich4-4  has  shown  that  the  pulling  effect  arises  from  a  summation  over  terms  with  the  coefficients 

expl/(/C^  +  Kl  -  K - -  K}):/2k}.  (23) 

from  which  the  Bragg  contributions  m  -  n,  m' —  n'  and  m  —  n',  m' —  n  are  excluded.  Here, 
K„  —  mirl  W  (where  m  is  an  integer)  is  the  eigenvalue  of  the  m th  waveguide  mode.  A  large  degree  of 
cancellation  will  occur  among  these  terms  if 

*-'(tfmax-  Kmm)  IKlz  >7T,  (24) 

where  AK  -  ir/W  is  the  adjacent  mode  spacing,  and  Ac  =  l/.?7is  the  distance  required  for  an  appreci¬ 
able  secular  change  in  the  backscatter  amplitude.  Noting  that  the  transverse  K  vector  spread  can  be 
expressed  as  -  ATmm  =  k0D,  one  finally  obtains  the  criterion 

9D/g7  >  tV  (25) 

for  the  pulling  effect  to  be  small.  In  the  simulations  performed  here,  W  —  2  mm;  hence  (25)  requires 
k9D/gl  >  kW  —  1.8  x  104,  in  good  agreement  wuh  Figs.  7  and  9.  Condition  (25)  can  be  rewritten  in 
a  useful  form  by  including  the  gain  requirement  gIL  >  10  to  obtain 

0£  >  (glL)  W/L  >  10W/L.  (26) 

The  second  notable  result  in  Sec.  Ill  confirms  that  the  effective  beam  divergence  0O  is  a  combina¬ 
tion  of  the  aberrator  divergence  9A  and  the  contribution  from  phase  curvature  introduced  by  a  lens.  If 
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a  large  initial  phase  curvature  is  present,  one  can  therefore  obtain  good  conjugation  fidelity  in  an  optical 
waveguide  with  weakly  aberrated  (or  even  unaberrated)  beams,  as  long  as  the  incident  profile  fills  the 
aperture.  It  is  worth  noting  here  that  the  incident  phase  curvature  also  plays  an  important  role  in  the 
focused  configuration.  For  example,  a  beam  of  width  WA  traversing  an  aberrator  of  divergence  #  ,  and 
contiguous  lens  of  focal  length  / will  produce  small  scale  intensity  inhomogeneities  characterized  by  ray 
crossing  angles  0  <  0F  -  WJf  near  the  focal  region.  Since  90  =  9f,  one  can  usually  satisfy  the  con¬ 
dition  kO'p  »  ,?/  (which  should  be  applicable  to  continuum  modes)  with  only  moderately  aberrated 
beams  as  long  as  the  lens  F  number  is  kept  low.  The  important  role  played  by  phase  curvature  has  also 
been  noted  in  earlier  work,’5  where  the  incident  beam  modulation  was  entirely  ir  the  amplitude  rather 
than  phase. 

In  Sec.  IV,  the  conjugation  fidelity  was  found  to  improve  with  pump  depletion  in  all  of  the  cases 
considered.  This  result  is  consistent  with  the  earlier  simulations  in  the  focused  configuration.11 
although  the  changes  are  less  dramatic  because  the  initial  fidelities  are  generally  higher. 

Although  fidelities  as  high  as  98%  were  found,  the  near  field  inhomogeneities  (^2:1  peak  to  val¬ 
ley)  persisted  in  all  cases.  This  behavior  can  be  adequately  explained  by  a  simple  model  introduced  in 
Ref.  (1 1).  One  approximates  the  backscatter  amplitude  at  the  right  hand  side  of  the  aberrator  by  a  con¬ 
stant,  plus  contributions  from  a  limited  number  of  modes  with  random  phase  amplitudes  a  v;  i.e., 

£*(*)  =  R[/ 2  E*l\\  +  £  “v  exp(/27r.Vjt/  WA)\,  (27) 

l  .vTo  I 

where  £ t  is  constant  within  interval  —  WJ2  <  x  <  +  WJ2,  and  zero  outside.  This  gives  the  near 
field  intensity 

\Es(x)\2~  Rl£il2[l  +  £  <*v  SKpU2nNx/Wt)  +  c.c.j  (28) 

l  v*0  j 

in  the  case  of  interest  where  la  vl2  <  <  1,  and  the  fidelity  (from  (9)) 

H  =  [l  +  £  l«  vP  (29) 

,v~> 

For  Nu  occupied  modes  of  comparable  RMS  amplitude  aRMs,  these  expressions  yield  the  approximate 
results 

l£slix/l£sU,n  ®  <1  +  2/VV/2  aRMS)/n  -  2/V,'/2  «RMS>  <30> 

//  =  I  —  /Vv/  aRMs-  (31) 

Considering  the  example  of  Fig.  11,  we  have  ,\u  aRMs  ~  1  ~  H  -  .03,  giving  l£sUax/l^slmm  ~  2  in 
reasonable  agreement  with  the  near  field  structure  shown  in  the  figure.  It  is  interesting  to  note  that 
even  for  H  -*  99%,  this  model  predicts  near  field  intensity  variations  up  to  1.5:1. 
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Fig.  2  —  Simulation  of  SBS  wavefront  conjugation  in  an  optical  waveguide  with  a  75  cm  gain  length,  showing  intensity  protiles  of 
the  incident  light  (solid  lines)  and  backscatter  (dashed  lines)  with  negligible  pump  depletion.  (The  incident  intensity  is  uniform 
across  the  aberrator.) 
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Fig.  4  —  Conjugation  fidelity  H and  conjugate  mode  enhancement  factor  fop  lEqs.  116)  and  <  I9)| 
vs  average  (amplitude)  gain  exponent  for  the  same  uberrator  and  gain  length  used  in  Fig.  2 
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Fig.  6  —  Conjugation  fidelity  vs  kSfrl gl.  as  calculated  from  Ref.  5  (dashed  curve)  and  simulations  with  aberrators  of  10  mrad 
(solid  diamonds)  and  20  mrad  (open  diamonds)  far  field  divergence.  Each  point  shows  the  average  and  mean  deviation  of  live 
simulations  using  independent  aberrators. 
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Fig.  10  —  Conjugation  fidelity  H  and  reflectivity  R  vs  incident  intensity  from  simulations  of  the  model  shown  in  Fig.  8  (L  -  75 
cm,  ■*  5  mrad,  0F  “  20  mrid).  The  dashed  line  shows  the  fidelity  in  the  absence  of  pump  depletion. 
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